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The present invention relates to a method and apparatus 
for analysing movement of a golf ball. 

Systems are known for tracking golf balls hit by a golf 
club. For example, US patent No. 617972 0 discloses an 
arrangement that matches the impact location of an object 
at or near a predetermined target with one of many shots 
at the target. US patent No. 54 01026 relates to a 
Doppler radar system which is used to calculate the carry 
distance of a golf shot. 

Radar devices which utilise the Doppler frequency shift 
to measure the velocity of a moving object are known in 
many applications. With such devices electromagnetic 
energy, such as microwave radar energy, is transmitted 
towards and reflected from a moving object. The energy 
undergoes a frequency shift, the magnitude of which is 
proportional to the velocity of the object relative to 
the transmitter. Samples of the transmitted and 
reflected radiation are mixed and processed to obtain a 
difference signal having a frequency which is equal to 
the difference between the transmitted and reflected 
frequencies, this difference being the Doppler shift. 
Once the difference frequency has been obtained, the 
relative velocity of the object can be readily calculated 
using known methods. 

In addition to using Doppler radar for measuring 
characteristics -of golf ball travel, other techniques 
have been devised, such as Doppler laser, sonar, light 



2 



gates and cameras, all having various degrees of 
accuracy. The most common approach is to use Doppler 
radar to measure the apparent speed of the golf ball 
directly from or towards the Doppler radar device. 
5 However, known Doppler radar devices suffer from a major 
problem in that no allowance is made for any deviation 
from the axis extending from or towards the radar device . 
This measurement of speed of the golf ball directly from 
or towards the Doppler radar device is hereinafter 
10 referred to as the "apparent velocity" of the golf ball. 
The measurement of the apparent velocity of a golf ball 
at launch, i.e. immediately after striking by the golf C ■ 

club, will provide only a very crude estimation of how 
far the golf ball will travel (the "carry distance"). 

It is an object of the present invention to provide a 
device with improved golf ball tracking. 

According to a first aspect of the present invention, 

2 0 there is provided a method of analysing movement of a 

golf ball, the method including transmitting radiation 
towards a location expected to lie in the path of a golf 
ball; receiving reflected radiation and generating a ( ) 

signal representative thereof; processing the signal to 
25 derive a value indicative of the velocity of the golf 

ball in a predetermined direction; processing the signal 
to derive a value indicative of the angular deviation of 
the golf ball trajectory from the predetermined 
direction; and deriving from the velocity and deviation 

3 0 values an indication of the speed of the ball along its 

trajectory. 
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In the embodiment, the derivation of a value indicative 
of the velocity of the golf ball in a predetermined 
direction can be considered to correspond to the 
"apparent velocity" mentioned above. By additionally 
5 processing the signal to derive a value indicative of 
angular deviation of the golf ball trajectory from the 
predetermined direction as well, a much better indication 
of the speed of the golf ball along its trajectory can be 
derived. 

10 

In the embodiment the velocity of the golf ball in the 
predetermined direction is determined by measuring the 
Doppler shift in the radiation. 

15 In the embodiment the step of receiving reflected 

radiation includes receiving radiation from at least two 
directions parallel to an axis. The step of receiving 
radiation advantageously further includes receiving 
radiation from at least one direction parallel to a 
2 0 further axis, different from the other axis. The two axes 
may be parallel. 

The method advantageously includes analysing a signal 
generated from the received radiation and determining the 
angular deviation of the golf ball from the or each axis. 
By determining the angular deviation of the golf ball 
from two different axes, the speed and direction of golf 
ball travel can be calculated relatively accurately. 
Determining to angular deviation may include determining 
the path length difference of the received radiation, and 
prefer ably- -include s -determining— the -phase- di-f f-er-enoe of 
the received radiation. 



25 



30 
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Means for receiving radiation may be provided, being 
sensitive to radiation substantially within a 20° 
vertical field and a 10° horizontal field. The means for 
5 receiving radiation may be positioned substantially 190mm 
above ground level and between 3.3 and 3 . 5m. behind the 
expected golf ball launch position. 

The method preferably also includes detecting when the 
10 golf ball is struck by a golf club. 

O 

According to the embodiment, using accurate angular 
measurements along the ball trajectory, it is possible to 
correct the apparent speed of the golf ball obtained from 

15 the Doppler shift measurement to true speeds at any point 
along the trajectory. Regardless of the club used, or the 
prevailing atmospheric conditions, the actual trajectory 
is either measured in full, or the device characterises 
the dynamic behaviour of a particular golf ball during 

2 0 flight and hence making it possible to predict to within 
a couple of metres, reliably the carry distance of a golf 
ball. 

O 

According to a second aspect of the present invention, 

2 5 there is provided apparatus for analysing movement of a 

golf ball, the apparatus including means for transmitting 
radiation towards a location expected to lie in the path 
of a golf ball; means for receiving reflected radiation 
and generating a signal representative thereof; means for 

3 0 processing the signal to derive a value indicative of the 

velocity of the goTf "ban ln"a-^ — 

means for processing the signal to derive a value 
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indicative of- angular deviation of the golf ball 
trajectory from the predetermined direction; and means of 
deriving from the velocity and deviation values an 
indication of the speed of the golf ball along its 
5 trajectory. 

For a better understanding of the present invention, an 
embodiment will now be described by way of example, with 
reference to the accompanying drawings, in which 

10 

Figure 1 shows a perspective view of a golf radar 
enclosure according to the present invention; 

Figure 2 shows a front elevational view of the golf radar 
15 enclosure of Figure. 1; 

Figure 3 shows schematically how the golf radar enclosure 
will be positioned with respect to a golfer as he strikes 
the ball; 

20 

Figure 4 shows schematically the relative positions of 
the transmitter and receiver elements of the system; 

Figure 5 shows an electronic functional block diagram of 

2 5 the system; 

Figure 6 illustrates the principle for measuring angular 
deviation of the golf ball; 

3 0 Figures 7 and 8 show process flow diagrams for tracking 
rhe^ golf "ball y ~ - - 
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Figure 9 shows a plot of an example signal received by a 
receiver of the system against time; 

Figure 10 shows a plot of the velocity waterfall signal 
5 derived from one of the receivers; 

Figure 11 shows a plot of the measured trajectory of the 
golf ball in elevation and azimuth against time; 

10 Figure 12 shows a plot of the velocity waterfall signal 
indicating golf club movement; and 

Figure 13 shows a golf ball trajectory. 

15 In the figures like elements are generally designated 
with the same reference numeral. 

Figures 1 and 2 show the golf radar enclosure 1. The 
enclosure l comprises front 3 and rear 5 main surfaces. 

2 0 The front surface includes a display unit 7, such as a 

liquid crystal display unit for displaying calculated 
data concerning a particular shot. A transmitter 9 and 
three receivers 11, 13 and 15 are positioned on the front 
surface 3'. The enclosure l includes digital signal 
25 processing electronics, display processing electronics 

and a power supply (none of which are shown in Figures 1 
and 2) . 

As shown in Figure 3, the enclosure 1 is positioned 

3 0 approximately 190mm above ground level and approximately 

3T3~Xo~ 3T5~"me t res BeKlnS ' "^S^osiT:"ion" of ' a ~ goTf er "17 ~ who ~ 
is about to strike a golf ball, for example as he or she 
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tees off. The enclosure 1 is positioned directly behind 
the golf ball launch position. The enclosure 1 is 
oriented so that the front surface 3 faces the golf ball 
launch position. 

Figure 4 shows schematically the arrangement of the 
transmitter 9 and the three receivers 11, 13 and 15 of 
the enclosure 1 . Energy is transmitted towards the golf 
ball launch position by transmitter 9. The energy is 
reflected from every object within the field of view of 
the transmitting device 9. The golfer 17 striking a ball 
together with his or her club 19, as well as the ball 
struck and balls from other players are all potential 
reflectors of the transmitted energy. The receivers 11, 
13 and 15 are synchronous receivers and measure the 
reflected energy from all targets within their fields of 
view. The transmitter and three receivers are generally 
placed in a particular pattern to optimise the signal 
quality from the golf ball and to minimise the negative 
influence from the golfer 17, swinging club 19 and any 
other sources of reflection. 

The receivers 11, 13 and 15 are spaced apart by a 
predetermined distance. In order to measure the angular 
deviation of the golf ball from the vertical, receivers 
13 and 15 are spaced apart vertically by a distance dv. 
In order to measure the deviation from the horizontal, 
the receivers 11 and 13 are spaced apart horizontally by 
a distance dn. The receiver pair 11, 13 have a sensor 
axis 2 0 extending out of the page in Figure 4, and the 

receiver—pai-r- 13-, — 3.-5— have— a— sensor— axis— 2-i-—- al so 

extending out of the page in Figure 4 . 
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The antenna beamwidths are approximately 2 0 degrees in 
vertical and 10 degrees in the horizontal direction. The 
antenna beam widths used, as well as the placement 
relative to the golfer and the ground, are selected in 
order to minimise the effect of multipath and hence 
distortion of the relative phase measurements, which in 
turn reflect in the angular measurements. 

It will be noted that, in this embodiment, the second 
receiver 13 contributes to angular deviation data from 
the horizontal and the vertical. This is a cost-effective 
arrangement. However, the second receiver 13 could be 
replaced with two separate receivers, a first receiver 
being dedicated to measure horizontal deviation with the 
receiver 11, and a second receiver dedicated to measuring 
vertical deviation with receiver 15. 

Figure 5 shows schematically the components contained 
within the enclosure 1. A microwave oscillator 22 
generates a signal of predetermined frequency, which is 
amplified by power amplifier 23. The amplified signal is 
passed to both the transmitter antenna 9 and to 
respective mixers 25, 27 and 29 by power splitter 31. As 
described above, reflected radiation is received by 
receiver antennae 11, 13 and 15. The signals from the 
receiving antennae 11, 13 and 15 are passed to respective 
low noise amplifiers 33, 35 and 27. The amplified 
signals are applied respectively to mixers 25, 2 7 and 29, 
where they are mixed with the amplified oscillating 

""signal:— f rour ^the-beam splirtter--31- — The -Ttixxed" -analogue 

signals are then converted to digital signals by analogue 



to digital converter 39 and subsequently processed by 
signal processor 41. The outputs from the digital signal 
processor 41 are then configured by display processor 4 3 
for presentation on the display 7 of the enclosure 1. 

5 

The operation of the apparatus will now be described. 
The description of how to obtain the "apparent velocity" 
will, however, be omitted for the sake of brevity, as 
this would be known to a person skilled in the art . As 
10 mentioned above, the apparent velocity differs from the 
actual speed. This difference can be described by the 
following formula: 

Velocity (apparent) = Velocity (actual) x 
15 cos (vertical ball angle - device pointing angle) x 

cos (horizontal ball angle - device pointing angle) 

For example, for a ball launched at 3 0° vertical deviation 
from axis 21 and 20° in azimuth deviation from axis 20 
2 0 (the device pointing angle being 0° - zero in both cases) , 
the error in the apparent velocity compared to the actual 
velocity would amount to approximately 18.7%. 

The vertical ball angle is derived from the signals 
2 5 received at the main receiving antenna 13 and the 

vertical receiving antenna 15. The horizontal ball angle 

is derived from the signals received by the main 

receiving antenna 13 and horizontal receiving antenna 11. 

The interf erometry principle used to track the ball will 
30 be described with reference to Figure 6, and is explained 

by the publication "Introduction to Monopulse" , Donald R. 

Rhodes, PhD, which is incorporated herein by reference. 
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Figure 6 shows two receivers, which could correspond to 
receiver pair 11, 13 or receiver pair 11, 15. The 
receivers are spaced apart by a distance D and are 
directed along an axis 45, which correspond to axis 2 0 or 
21. When the golf ball 4 7 is in a position deviated from 
the axis 45, the radiation Rl and R2 reflected from the 
golf ball 4 7 and received by "Receiver 1" and "Receiver 
2 M respectively will travel different distances to each 
of the respective receivers, the difference being AR. 

Energy wave path length difference = 2 AR = D.sinG 

The path length difference can be expressed as a 
proportion of the electrical wavelength 

Acp = (2nD/X) .sinG 

where Acp is the phase difference between the 
receiving channels, 

?l is the wavelength of the energy transmitted 
by the transmitter 9 to determine the angles to 
the ball 47 trajectory, 

0 is the angle in space relating to the 
electrical phase difference Acp. 

At any particular point the path length difference is a 
direct indication of the angle of arrival of the plane 
wave front, indicating the angle of deviation e from the 
axis - 45- -of -the -receivers-- 



11 

Figures 7 and 8 show the processes performed by the 
digital signal processor 41 and display processor 43. 
When the digital signal processor 41 is initiated it 
synchronously samples received data from each of the 
5 receivers 11, 13 and 15 at step 50. Figure 9 shows a 

typical time signal series as measured by any one of the 
three receivers 11, 13 or 15. By sampling the time 
signals from each of the receivers 11, 13 and 15 
synchronously the relative phase information is preserved 

10 .in the signals. Data is collected and recorded from the 
receivers 11, 13 and 15 continuously and simultaneous 
digital signal processing is performed on the recorded 
data to provide velocity and angular information of items 
within the field of view of the sensors at least every 

15 two to 15ms . 

A 64 point fast furrier transform (FFT) is performed on 
the signals from the three receivers 11, 13 and 15 in 
real time to produce a velocity spectrum of all moving 

2 0 objects within the field of view of the receivers (step 

52) . At this stage the golf ball has not been struck and 
the digital signal processor 41 analyses the velocity 
spectra to identify when a golf club swing occurs, which 
is the trigger for measuring golf ball movement. The 

2 5 magnitude and tempo of a typical golf club swing is 
characteristic for a particular set up, which is 
mathematically modelled and used as a unique 
identification in order to eliminate triggering from 
neighbouring club swings, as well as balls which may come 

30 within the field of view of the receivers 11, 13 and 15. 

The^magnilTudes - of threat 'ime~~s ignal"of — alrl— t hr e e -rec e ivers — 
are compared with a calibrated range of values to ensure 
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a valid trigger condition. In addition, the apparent 
velocity profile of the club is used a characteristic to 
validate the trigger condition. Finally, the angles both 
vertical and horizontal of the club swing are analysed to 
further validate the swing path as a conclusive method of 
verifying that a club swing has in fact occurred. Figure 
10 shows typical club strike speed behaviour in the 
format of a velocity waterfall series. The trigger 
mechanism relies on the fact that the club's apparent 
rate as determined by the digital signal processor 41 is 
characteristic and lies between a minimum and maximum 
rate. The signal processor 41 measures the acceleration 
and deceleration of the club as it moves through the 
swing plane. If a golf ball is struck, the club as well 
as the ball appears as two moving objects, with the ball 
speed generally much higher than the club speed. The 
angular movement of the club face within the antenna 
beams is also measured, and hence a confirmation of the 
club swing trajectory is readily made as confirmation 
that a club swing has indeed taken place, right in front 
of the radar of concern and not introduced by 
neighbouring golfers, etc. 

If the digital signal processor 41 determines that a club 
swing has not occurred at step 54, the digital signal 
processor 41 returns to step 5 0 and monitor for a club 
swing. However, if it is determined that a club swing 
has occurred the digital signal processor 41 performs 
step 56 which analyses velocity waterfall signals 
computed for each of the three receivers 11, 13 and 15. 
The "^ela^^e^ ea cH~ "of T5e"tSree 

receivers is exactly the same. However, the phase 
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difference is extracted in the frequency domain and used 
for angular calculations (described later) . Figure 10 
shows a typical velocity waterfall plot. If the analysis 
of the velocity waterfall indicates that ball movement is 
5 not present the step 58 returns processing to step 50, 
where the signals from the receivers 11, 13 and 15 are 
sampled. However, if it is detected that ball movement 
has occurred step 60 identifies the start of the golf 
ball trajectory to. Step 62 then performs a 2 56 point FFT 

10 on the signals from the receivers 11, 13 and 15. An 

algorithm is used to filter out values calculated not to 
belong to the golf ball trajectory. At step 64 a 1024 
point FFT is performed to extract the phase difference 
between the receiving channels 11, 13 and 15. At step 66 

15 the deviation from the axis 21 of the receiver pair 13 

and 15 is calculated. At step 68 the deviation from the 
axis 2 0 of the receiver pair 11, 13 is calculated. The 
calculations at step 66 and 6 8 are performed using the_ 
principle described in relation to Figure 6 various 

20. distances from the enclosure 1. Doppler shift 

information is also extracted in the conventional manner, 
and at step 7 0 the apparent velocity obtained by the 
Doppler shift method is used to provide actual velocity 
data using the angular trajectory information obtained at 

25 step 66 and 68. 

Referring now to Figure 8, a mathematical model is used 
to extract the dynamic drag coefficient Cdo at step 72 and 
the dynamic lift coefficient CLo at step 74. 

30 

— : — The— dynamic —drag -and— lift— coeffix; i-ent-s— embody— fefae—ba-3ri 

physical properties, the wind and other atmospheric 



BNSDOCID: <GB 2380682A_I_> 



14 

conditions. The derivation of these coefficients allows, 
regardless of the club used and the prevailing 
atmospheric conditions, the dynamic behaviour of a 
particular golf ball strike to be characterised, hence 
making it. possible to predict reliably to within a couple 
of metres the carry distance of a golf ball. 

Ideally, the receivers 11, 13 and 15 would receive data 
relating to the entire golf ball trajectory from the 
initial strike by the club until the golf ball travel is 
finished. However, in practice, it is often only 
possible to track part of the trajectory. At step 76 the 
data obtained thus far is compared to an imperical model 
to compute the full trajectory of the golf ball when data 
is only available relating to part of the trajectory. 

The actual distance that a golf ball travels through air 
is related, inter alia, to the launch speed in the 
direction of movement, the vertical and horizontal launch 
angles, the spin rate and direction (spin vector) , the 
mass of the ball, the dynamic drag coefficient Cdo and the 
dynamic lift coefficient CLo (taking into account the 
current actual spin vector, wind, air temperature and 
other atmospheric conditions) . 

At step 78 the carry distance of the ball is displayed on 
the display 7 of the enclosure 1. At step 8 0 the launch 
velocity of the ball is displayed. At step 82 the club 
strike speed is displayed 

This* "data "is "displayed" u^uaTXy ^ 
ball being struck. 
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The above processes will now be further described with 
reference to Figure 13 . 

5 Step 1 

Process raw data to determine radial speed,. Vr, and 
elevation angle, <|>, measured by the radar. 

Step 2 

10 Use two first order polynomial fittings on Vr to find the 
^ starting point and to extract the good data (i.e. remove 

values not belonging to the golf ball trajectory) . 

Step 3 

15 Extract the real speed, V, and the angle, 9, by using an 
iteration formula. To initiate the process, use the 
initial values 9 = 18° and v = Vr and then iterate 3 0 
times : 

(t = time elapsed since to, g is gravity) 

20 

o 
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Ar = vA: 

Ax = Ar cos(&) 

Ay = Arsm(&) 







y - A 




6 = axctan 


A x 


r + ^ 




A x 


( y - k] 




I X + * J J 






v = v p sec(0 


-<0 





Step 4 

Extract the drag coefficient, and the lift 
coefficient, y, by fitting 5" order polynomials to the 
horizontal and vertical speed components. The speed 
components are given by x =vcos(0) and y = vsin(0) - 

Once the polynomial coefficients {*0 and {&} such that 
ir(r) = a, +2cr 2 /-r3a 3 r : -■-... + 6c:,,/ 5 

are known, calculate the coefficients with the formula 
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A 1 








/V 2 





o 
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_3.0_ 



Step 5 

Estimate horizontal distance travelled. This function 
uses a numerical solution of the model. 

x = -jjx* -yy 

y = -My\S>\- 

with the difference formulas 

At 



y*- g 



x M - 2x, + = -p(x,^ - x, \x, - r ; _, )-y - J>,_, ) 



which can be written in matrix form as 
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O 



10 



15 



20 



A 



where 



2 



b 2 = 2 v, - ,y M +My,\y i - | — — x M - 
The initial values v and 0 are needed to estimate xi and 
yi as follows : 

x, s=a,A/ + a 2 Ar 2 



a, = vcos(#) 



where 

Figure 11 shows the vertical and horizontal trajectory 
information as a function of time. This trajectory 
information together with a velocity waterfall of Figure 
10 is used to compute the ball trajectory and all 
associated ballistic perimeters, required to model the 
ball movement through space. 



2 5 The device of the embodiment uses microwave frequency 
radiation, but any electromagnetic, acoustic or other 
wave -like energy could be used. 



BNSDOCID: <GB_ 



_2380682A_I_> 



18 

CLAIMS: 

1. A method of analysing movement of a golf ball, the 
method including transmitting radiation towards a 
5 location expected to lie in the path of a golf ball; 
receiving reflected radiation and generating a signal 
representative thereof; processing the signal to derive a 
value indicative of the velocity of the golf ball in a 
predetermined direction; processing the signal to derive 
10 a value indicative of angular deviation of the golf ball 
trajectory from the predetermined direction; and deriving 
from the velocity and deviation values an indication of 
the speed of the golf ball along its trajectory. 

15 2. A method according to claim 1, wherein the velocity 
of the golf ball in a predetermined direction is 
determined by measuring Doppler shift of the radiation. 

3. A method according to claim 1 or 2 , wherein the step 

2 0 of receiving radiation includes receiving radiation from 

at least two directions parallel to an axis. 

4 . A method according to claim 3 , including receiving 
radiation from at least one direction parallel to a 
25 further axis, different from the other axis. 

5. A method according to claim 3 or 4 , including 
analysing a signal generated by the received radiation 
and determining the angular deviation of the golf ball 

3 0 from the or each axis. 
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6. A method according to claim 5, including determining 
the path length difference of the received radiation. 



C) 



( J 



7. A method according to claim 5 or 6, including 
5 determining the phase difference of the received 

radiation. 

8. A method according to any one of claims 1 to 7 , 
wherein means for receiving radiation is provided, being 

10 sensitive to radiation substantially within a 2 0° 
vertical field and a 10° horizontal field. 

9. A method according to any one of claims 1 to 8, 
wherein means for receiving radiation is provided, being 

15 substantially 190mm above ground level and between 3.3 
and 3,5m behind the expected golf ball launch position. 

10. A method according to any one of claims 1 to 9, 
including detecting when the golf ball is struck by a 

2 0 golf club. 

11. A method of analysing movement of a golf ball, 
substantially as hereinbefore described with reference to 
and/or as illustrated in any one of or any combination of 

2 5 the accompanying drawings. 

12. Apparatus for analysing movement of a golf ball, the 
apparatus including means for transmitting radiation . 
towards a location expected to lie in the path of a golf 

30 ball; means for receiving reflected radiation and 

generating -a -signal- represented ive - thereof ->•— me^ans-f- or- 

processing the signal to derive, a value indicative of the 
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velocity of the golf ball in a predetermined direction; 
means for processing the signal to derive a value 
indicative of angular deviation of the golf ball 
trajectory from the predetermined direction; and means of 
deriving from the velocity and deviation values an 
indication of the speed of the golf ball along its 
trajectory. 

13. Apparatus according to claim 12, including means for 
measuring Doppler shift of the radiation, to determine 
the velocity of the golf ball in a predetermined 
direction. 

14. Apparatus according to claim 12 or 13, wherein the 
receiving means receives radiation from at least two 
directions parallel to an axis. 

15. Apparatus according to claim 14, wherein the 
receiving means receives radiation from at least one 
direction parallel to a further axis, different from the 
other axis . 

16. Apparatus according to claim 14 or 15, including 
means for analysing a signal generated by the received 
radiation and determining the angular deviation of the 
golf ball from the or each axis. 

17. Apparatus according to claim 16, including means for 
determining the path length difference of the received 
radiation . 



18. Apparatus according to claim 16 or 17, including 
means for determining the phase difference of the 
received radiation. 

19. Apparatus according to any one of claims 12 to 18 , 
wherein means for receiving radiation is provided, being 
sensitive to radiation substantially within a 20° 
vertical field and a 10° horizontal field. 

20. Apparatus according to any one of claims 12 to 19, 
wherein means for receiving radiation is provided, being 
substantially 190mm above ground level and between 3.3 
and 3.5m behind the expected golf ball launch position. 

21. Apparatus according to any one of claims 12 to 20, 
including means for detecting when the golf ball is 
struck by a golf club. 

22. Apparatus for of analysing movement of a golf ball, 
substantially as hereinbefore described with reference to 
and/or as illustrated in any one of or any combination of 
the accompanying drawings. 
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